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Technology Focus: Test & Measurement 


Wireless Measurement of Contact and Motion Between 
Contact Surfaces 

A magnetic-field-response contact sensor is used. 

Langley Research Center, Hampton, Virginia 


This method uses a magnetic-field- re- 
sponse contact sensor that is designed to 
identify surface contact and motion be- 
tween contact locations. The sensor has 
three components: (1) a capacitor-induc- 
tor circuit with two sets of electrical con- 
tact pads, (2) a capacitor with a set of 
electrical contact pads, and (3) an induc- 
tor with a set of electrical contact pads. A 
unique feature of this sensor is that it is 
inherently multifunctional. Information 
can be derived from analyzing such sen- 
sor response attributes as amplitude, fre- 
quency, and bandwidth. A change in one 
attribute can be due to a change in a 


physical property of a system. A change 
in another attribute can be due to an- 
other physical property, which has no re- 
lationship to the first one. 

The sensor is powered and interro- 
gated without physical connection to a 
power source, microprocessor, data acqui- 
sition equipment, or electrical circuitry. It 
works with the magnetic-field-response 
recorder described in “Magnetic-Field-Re- 
sponse Measurement-Acquisition Sys- 
tem,” NASA Tech Briefs, Vol. 30, No. 6 
(June 2006) , page 28. 

The sensor (capacitor-inductor cir- 
cuit) is placed on the moving object as 



Edge View of Object is shown moving away from surface A toward surface B. Surface A has a capaci- 
tor with electrical contact pads. Surface B has an inductor with contact pads. The moving object has a 
inductor-capacitor circuit with two sets of contacts pads to complete the electrical connection with ei- 
ther surface. 


shown in the figure. When contact is 
made between stationary surface A, the 
capacitor and the capacitor-inductor cir- 
cuit form a circuit the response fre- 
quency of which is 

1 

oj = 

The contact pads could be compressible 
or spring-loaded electrical contacts that 
are electrically connected to the sensor. 
Movement away from surface A to the 
next surface results in the sensor re- 
sponse frequency shifting to 
1 

uj — 

The response amplitude decreases as 
the sensor moves away from the an- 
tenna. Contact with surface B shifts the 
sensor response frequency to 

__L HT 

2k \ LC 

Applications for this include being able 
to tell whether doors or hatches are 
sealed, surface bonds are secure (tile 
bonds, rubber bond to steal belts of tires, 
etc.), general knowledge of contact be- 
tween two surfaces, or separation of sur- 
faces. 

This work was done by Stanley E. 
Woodard of Langley Research Center and 
Bryant D. Taylor of Swales Aerospace. Fur- 
ther information is contained in a TSP (see 
page 1). LAR- 16849-1 




Wireless Measurement of Rotation and Displacement Rate 


A magnetic field response ser 

Langley Research Center, Hampton, 

A magnetic field response sensor is 
designed to measure displacement or 
rotation rate without a physical connec- 
tion to a power source, microprocessor, 
data acquisition equipment, or electri- 
cal circuitry. The sensor works with the 


is used in these measurements. 

inia 

magnetic field response recorder, which 
was described in “Magnetic-Field-Re- 
sponse Measurement-Acquisition Sys- 
tem,” NASA Tech Briefs, Vol. 30, No. 6 
(June 2006), page 28. These sensors are 
wirelessly powered and interrogated, 


and the measurement acquisition sys- 
tem and sensors are extremely light- 
weight. 

The response recorder uses oscillating 
magnetic fields to power the sensors. 
Once powered, the sensors respond with 


NASA Tech Briefs, November 2007 


5 





Two Response Sensors with different resonant 
frequencies are sown inside of wheel. 


their own magnetic field. For displace- 
ment/rotation measurements, the re- 
sponse recorder uses the sensor’s re- 
sponse amplitude, which is dependent 
on the distance from the antenna. The 


recorder’s antenna orientation and posi- 
tion are kept fixed, and the sampling pe- 
riod is constant. 

A sensor with fixed frequency and 
fixed orientation with respect to the re- 
sponse recorder antenna can be used 
for position and displacement measure- 
ments. If the sensor’s orientation is not 
fixed, but its trajectory is known, it may 
be possible to calibrate the response 
amplitude variation with trajectory. For 
rotational motion such as wheel speed, 
identifying the number of times the re- 
sponse amplitude exceeds threshold 
amplitude in a fixed time duration can 
be used to determine rotation rate. A 
wheel speed sensor is shown in the fig- 
ure. The sensor is a thin-film circuit 


placed inside the wall of a tire. As the 
sensor approaches the antenna, the am- 
plitude increases. The amplitude peaks 
at the closest point to the antenna then 
decreases producing changes in ampli- 
tude that are cyclical. When two sen- 
sors with different respond frequencies 
are used inside the wheel, rotation di- 
rection can be determined by identify- 
ing which sensor’s amplitude increases 
first. In addition, there is no mechani- 
cal wear because no gears are used in 
the design. 

This work was done by Stanley E. Woodard 
of Langley Research Center and Bryant D. 
Taylor of Swales Aerospace. Further informa- 
tion is contained in a TSP (seepage 1 ). LAR- 
16848-1 


Portable Microleak-Detection System 

Heating or cooling of a vacuum seal enables testing over a wide temperature range. 

Langley Research Center, Hampton, Virginia 


The figure schematically depicts a 
portable microleak-detection system that 
has been built especially for use in testing 
hydrogen tanks made of polymer-matrix 
composite materials. (As used here, “mi- 
croleak” signifies a leak that is too small to 
be detectable by the simple soap-bubble 
technique.) The system can also be used 
to test for microleaks in tanks that are 
made of other materials and that contain 
gases other than hydrogen. Results of cal- 
ibration tests have shown that measure- 
ment errors are less than 10 percent for 
leak rates ranging from 0.3 to 200 
cmVmin. 

Like some other microleak-detection 
systems, this system includes a vacuum 
pump and associated plumbing for sam- 
pling the leaking gas, and a mass spec- 
trometer for analyzing the molecular 
constituents of the gas. The system in- 
cludes a flexible vacuum chamber that 
can be attached to the outer surface of a 
tank or other object of interest that is to 
be tested for leakage (hereafter de- 
noted, simply, the test object). The gas 
used in a test can be the gas or vapor 
(e.g., hydrogen in the original applica- 
tion) to be contained by the test object. 
Alternatively, following common prac- 
tice in leak testing, helium can be used 
as a test gas. In either case, the mass 
spectrometer can be used to verify that 
the gas measured by the system is the 
test gas rather than a different gas and, 
hence, that the leak is indeed from the 
test object. 


Liquid-Hydrogen or Liquid-Helium 
Manifold in Test Object 


Test Object 


Capillary Tube- 



Valves. 




/ 


Vacuum Membrane 


\ 

Seal Containing 
Heating/Cooling Tube 


Mass 

Spectrometer 


Pressure and 
Temperature Sensors 



Vacuum Pump 


Control Volume 


The Portable Microleak-Detection System includes components in common with prior microleak-detec- 
tion systems, plus a seal-heating/cooling subsystem that enables testing over a wide temperature 
range. 


The flexibility of the chamber makes 
it adaptable to test objects having a vari- 
ety of simple or complex shapes. The 
flexible vacuum chamber includes an 
aluminized polyethylene terephthalate 
vacuum membrane that is sealed to the 
outer surface of the test object by a flex- 
ible, adhesive seal material. A scrim is 
placed between the inner surface of the 
membrane and the outer surface of the 
test object to maintain a gap to accom- 
modate the flow of any leaking gas. A 
capillary tube that passes through the 


seal connects the gap volume with the 
plumbing that leads to the mass spec- 
trometer, the vacuum pump, and a con- 
trol volume described next. 

The control volume has a known size 
and is instrumented with pressure and 
temperature sensors. In use, the control 
volume is evacuated, then disconnected 
from the vacuum pump, and then the 
pressure and temperature are measured 
as the leaking gas flows into the control 
volume. By use of the ideal-gas law, the 
rate of leakage can be calculated from 
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the temperature and measured rate of 
increase of pressure. 

An unusual feature of this system is a 
heating/cooling subsystem that in- 
cludes a tube embedded in the flexible 
adhesive seal. A heating or cooling liq- 
uid can be circulated through this tube 
to maintain the seal at or near room 
temperature, where it is most effective, 


regardless of the temperature of the 
test object or the environment. The 
heating/cooling subsystem is essential, 
for example, for maintaining an effec- 
tive seal for testing a tank, pipe, valve, 
or other object that contains liquid hy- 
drogen or other cryogenic fluid. The 
heated/cooled seal enables testing at 
temperatures from -455 to +350 °F 


(about -271 to +177 °C), even in the 
presence of distortions caused by me- 
chanical and thermal loads applied to 
the test object. 

This work was done by H. Kevin Rivers and 
Joseph G Sikora of Langley Research Center and 
Sankara N. Sankaran of Lockheed Martin 
Space Operations. Further information is con- 
tained in a TSP (seepage 1 ). LAR-16139 


# Free-to-Roll Testing of Airplane Models in Wind Tunnels 

Causes of, and cures for, wing-drop/rock behavior can be evaluated. 

Langley Research Center, Hampton, Virginia 


A free-to-roll (FTR) test technique and 
test rig make it possible to evaluate both 
the transonic performance and the wing- 
drop/rock behavior of a high-strength 
airplane model in a single wind-tunnel 
entry. The free-to-roll test technique is a 
single degree-of-motion method in which 
the model is free to roll about the longi- 
tudinal axis. The rolling motion is ob- 
served, recorded, and analyzed to gain in- 
sight into wing-drop/rock behavior. 

Wing-drop/ rock is one of several 
phenomena symptomatic of abrupt 
wing stall. FTR testing was developed as 
part of the NASA/Navy Abrupt Wing 
Stall Program, which was established 
for the purposes of understanding and 
preventing significant unexpected and 
uncommanded (thus, highly undesir- 
able) lateral-directional motions associ- 
ated with wing-drop/rock, which have 
been observed mostly in fighter air- 
planes under high-subsonic and tran- 
sonic maneuvering conditions. Before 
FTR testing became available, wing- 
rock/drop behavior of high-perform- 
ance airplanes undergoing develop- 
ment was not recognized until flight 
testing. FTR testing is a reliable means 
of detecting, and evaluating design 
modifications for reducing or prevent- 
ing, very complex abrupt wing stall 
phenomena in a ground facility prior 
to flight testing. 

The FTR test rig was designed to re- 
place an older sting attachment butt, 
such that a model with its force bal- 
ance and support sting could freely ro- 
tate about the longitudinal axis. The 
rig (see figure) includes a rotary head 
supported in a stationary head with a 
forward spherical roller bearing and 
an aft needle bearing. Rotation is am- 
plified by a set of gears and measured 
by a shaft-angle resolver; the roll angle 
can be resolved to within 0.067° at a 


rotational speed up to 
l,000°/s. An assembly of elec- 
trically actuated brakes be- 
tween the rotary and stationary 
heads can be used to hold the 
model against a rolling torque 
at a commanded roll angle. 
When static testing is required, 
a locking bar is used to fix the 
rotating head rigidly to the sta- 
tionary head. Switching be- 
tween the static and FTR test 
modes takes only about 30 
minutes. The FTR test rig was 
originally mounted in a 16-ft 
(~4.0-m) transonic wind tun- 
nel, but could just as well be 
adapted to use in any large 
wind tunnel. 

In one series of tests on the 
FTR rig, static and dynamic 
characteristics of models of four 
different fighter airplanes were 
measured. Two of the models 
exhibited uncommanded lateral 
motions; the other two did not. 
A figure of merit was developed 
to discern the severity of lateral 
motions. Using this figure of 
merit, it was shown that the FTR 
test technique enabled identifi- 
cation of conditions under 
which the uncommanded lat- 
eral motions occurred. The 
wind-tunnel conditions thus 
identified were found to be cor- 
related with flight conditions 
under which the corresponding 
full-size airplanes exhibited un- 
commanded lateral motions. 

This work was done by Francis J. 
Capone, D. Bruce Owens, and 
Robert M. Hall of Langley Research 
Center. Further information is con- 
tained in a TSP (see page 1). LAR- 
17133-1 
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EXPLODED VIEW SHOWING PARTS OF TEST RIG 


The FTR Test Rig allows free rotation of the airplane model 
around the centerline. Optionally, brakes can hold the 
model at a commanded roll angle, or the rig can be locked 
against rotation for conventional static testing. 
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# Cryogenic Shrouds for Testing Thermal-Insulation Panels 

These shrouds enable maintenance of required thermal and mechanical conditions. 

Marshall Space Flight Center, Alabama 


Cryogenic shrouds have been de- 
signed and built for use in thermome- 
chanical testing of samples of thermal- 
insulation panels on cryogenic vessels. 
In the original application for which 
these shrouds were specifically designed, 
the samples are representative of the 
large-area thermal-insulation panels on 
the space-shuttle external tanks that 
hold liquid hydrogen and liquid oxygen, 
and the purpose of the testing is to 
demonstrate the ability of bonded layers 
in the panels to resist delamination 
under a combination of applied uniaxial 
mechanical loads and realistic opera- 
tional temperatures. Presumably, the 
shrouds and the tests performed by use 
of them could be modified to enable 
similar evaluation of thermomechanical 
properties of thermal-insulation panels 
for cryogenic vessels other than the ex- 
ternal tanks of the space shuttles. 

The shrouds are required to enable 
maintenance of required temperatures 
on the inner and outer surfaces of the 
thermal-insulation-panel samples, to en- 
able visual observation of the outer sur- 


faces of the samples, and not to intro- 
duce any measurable loads into the pan- 
els. For each panel sample, there are two 
shrouds: one to be mounted on the 
inner surface (the surface that would be 
in contact with a tank containing a cryo- 
genic liquid during normal use) and 
one to be mounted on the outer surface 
(the surface that would be exposed to 
ambient air or other warmer environ- 
ment during normal use). The shrouds 
for testing specimens of thermal-insula- 
tion-panels for the liquid-hydrogen tank 
are made largely of titanium; the 
shrouds for testing specimens of ther- 
mal-insulation-panels for the liquid-oxy- 
gen tank are made largely of an alu- 
minum-lithium alloy. 

The specific temperature require- 
ments are the following: The inner 
shroud must make it possible to main- 
tain a temperature of -321 °F (-196 °C) 
[the approximate temperature of liquid 
nitrogen] or -453 °F (about -269 °C) 
[the approximate temperature of liquid 
helium] on the inner face of the sample. 
The outer shroud must make it possible 


to maintain a temperature between -30 
and 0 °F (between about -34 and about 
-18 °C) on the outer surface of the sam- 
ple by blowing a cryogenic gas or missile- 
grade air along that surface. To enable 
viewing of the outer surface of the sam- 
ple during testing, the outer shroud in- 
cludes a window comprising two layers 
of poly (methyl methacrylate) with a gap 
between them to reduce fogging. 

To ensure that the shrouds do not in- 
troduce any measurable loads into a 
panel specimen, the shrouds are cush- 
ioned on the specimen by seals made of 
a fluoropolymer-membrane / fabric com- 
posite material and are held in place on 
the specimen by means of symmetrically 
placed clamps with poly(tetrafluoroeth- 
ylene) pads. Instrumentation ports for 
thermocouples and strain gauges used 
in the tests are incorporated into the 
shrouds. 

This work was done by Jeffrey Norris, 
Robert Carroll, and Charles Kirch of Lock- 
heed Martin Corp. for Marshall Space Flight 
Center. Further information is contained in a 
TSP (see page 1 ). MFS-32444-1 
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Electronics/Computers 


Optoelectronic System Measures Distances to Multiple 
Targets 

This system operates at an update rate >10 Hz. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An optoelectronic metrology appara- 
tus now at the laboratory-prototype 
stage of development is intended to re- 
peatedly determine distances of as 
much as several hundred meters, at sub- 
millimeter accuracy, to multiple targets 


in rapid succession. The underlying 
concept of optoelectronic apparatuses 
that can measure distances to targets is 
not new; such apparatuses are com- 
monly used in general surveying and 
machining. However, until now such ap- 


paratuses have been, variously, con- 
strained to (1) a single target or (2) 
multiple targets with a low update rate 
and a requirement for some a priori 
knowledge of target geometry. When 
fully developed, the present apparatus 



This Apparatus Includes N Laser Assemblies (of which only the first and N th are shown here), for measuring distances to N targets. The lasers are turned 
on, one at a time, to illuminate their targets for short intervals in rapid succession to obtain a high update rate. 
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would enable measurement of distances 
to more than 50 targets at an update 
rate >10 Hz, without a requirement for 
a priori knowledge of target geometry. 

The apparatus (see figure) includes a 
laser ranging unit (LRU) that includes 
an electronic camera (photo receiver), 
the field of view of which contains all 
relevant targets. Each target, mounted 
at a fiducial position on an object of in- 
terest, consists of a small lens at the out- 
put end of an optical fiber that extends 
from the object of interest back to the 
LRU. For each target and its optical 
fiber, there is a dedicated laser that is 
used to illuminate the target via the op- 
tical fiber. The targets are illuminated, 
one at a time, with laser light that is 
modulated at a frequency of 10.01 
MHz. The modulated laser light is emit- 
ted by the target, from where it returns 
to the camera (photodetector), where 
it is detected. 

Both the outgoing and incoming 
10.01-MHz laser signals are mixed with 


a 10-MHz local-oscillator to obtain 
beat notes at 10 kHz, and the differ- 
ence between the phases of the beat 
notes is measured by a phase meter. 
This phase difference serves as a meas- 
ure of the total length of the path trav- 
eled by light going out through the op- 
tical fiber and returning to the camera 
(photodetector) through free space. 

Because the portion of the path 
length inside the optical fiber is not or- 
dinarily known and can change with 
temperature, it is also necessary to meas- 
ure the phase difference associated with 
this portion and subtract it from the 
aforementioned overall phase differ- 
ence to obtain the phase difference pro- 
portional to only the free-space path 
length, which is the distance that one 
seeks to measure. Therefore, the appara- 
tus includes a photodiode and a circula- 
tor that enable measurement of the 
phase difference associated with propa- 
gation from the LRU inside the fiber to 
the target, reflection from the fiber end, 


and propagation back inside the fiber to 
the LRU. Because this phase difference 
represents twice the optical path length 
of the fiber, this phase difference is di- 
vided in two before subtraction from the 
aforementioned total-path-length phase 
difference. 

Radiation-induced changes in the 
photodetectors in this apparatus can af- 
fect the measurements. To enable cali- 
bration for the purpose of compensa- 
tion for these changes, the apparatus 
includes an additional target at a known 
short distance, located inside the cam- 
era. If the measured distance to this tar- 
get changes, then the change is applied 
to the other targets. 

This work was done by Carl Christian Liebe, 
Alexander Abramovici, Randall Bartman, 
Jacob Chapsky, John Schmalz, Keith Coste, Ed- 
ward Litty, Raymond Lam, and Sergei Jerebets 
of Caltech for NASA’s Jet Propulsion Labora- 
tory. Lurther information is contained in a 
TSP (seepage 1 ). 

NPO-43240 


©Tachometers Derived From a Brushless DC Motor 

Neither excitation of the motor nor mechanical brushes is necessary. 

Marshall Space Flight Center, Alabama 


The upper part of the figure illus- 
trates the major functional blocks of a 
direction-sensitive analog tachometer 
circuit based on the use of an unexcited 
two-phase brushless dc motor as a rota- 
tion transducer. The primary advan- 
tages of this circuit over many older 
tachometer circuits include the follow- 
ing: 

• Its output inherently varies linearly 
with the rate of rotation of the shaft. 

• Unlike some tachometer circuits that 
rely on differentiation of voltages with 
respect to time, this circuit relies on in- 
tegration, which results in signals that 
are less noisy. 

• There is no need for an additional 
shaft-angle sensor, nor is there any 
need to supply electrical excitation to 
a shaft-angle sensor. 

• There is no need for mechanical 
brushes (which tend to act as sources 
of electrical noise) . 

• The underlying concept and electrical 
design are relatively simple. 

This circuit processes the back-elec- 
tromagnetic force (back-emf) outputs 
of the two motor phases into a voltage 
directly proportional to the instanta- 
neous rate (sign • magnitude) of rota- 


Multiplier Circuits 



Integrator 

Circuits 


MAGNITUDE-AND-SIGN CIRCUIT 



MAGNITUDE (ONLY) CIRCUIT 


These Analog Tachometer Circuits perform straightforward operations on the back-emf outputs of a 
brushless dc motor to generate voltages proportional to the rate of rotation of the shaft. 
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tion of the shaft. The processing in this 
circuit effects a straightforward combi- 
nation of mathematical operations lead- 
ing to a final operation based on the 
well-known trigonometric identity (sin 
x) 2 + (cos x) 2 = 1 for any value of x. The 
principle of operation of this circuit is 
closely related to that of the tachometer 
circuit described in “Tachometer De- 
rived From Brushless Shaft-Angle Re- 
solver” (MFS-28845), NASA Tech Briefs , 
Vol. 19, No. 3 (March 1995), page 39. 
However, the present circuit is simpler 
in some respects because there is no 
need for sinusoidal excitation of shaft- 
angle-resolver windings. 

The two back-emf signals are k 6 sin 6 
for phase A and k 6c os 6 for phase B, 
where k is a constant that depends on the 
electromagnetic characteristics of the 
motor, 6 is the instantaneous shaft angle, 
and the overdot signifies differentiation 
with respect to time. Note that 6 is the 
quantity that one seeks to measure. 

Each back-emf signal is fed to one of 
two inputs of a multiplier circuit of 
gain k 2 dedicated to its respective 
phase. Each of these signals is also inte- 
grated with a suitable time constant 
and gain to obtain a voltage of ^sin 6 
for phase A and -^cos 6 for phase B 
(where k x is a constant that incorpo- 


rates the combined effects of the gain 
and the time constant). The output of 
the integrator for phase B is inverted to 
obtain a voltage ^cos 6. Each of these 
signals is fed to the other input termi- 
nal of the multiplier circuit for its re- 
spective phase. 

The multiplier circuit for phase A thus 
generates an output signal proportional 
to both of its inputs; namely F, 0(sin Q) 2 , 
where ^ In a similar manner, the 

multiplier circuit for phase B generates an 
output signal of h, 6(c os 6) 2 . These signals 
are fed to an adder circuit. By virtue of the 
identity (sin 6) 2 + (cos 6) 2 = 1, the output 
of the adder is simply k s 6. 

The lower part of the figure illustrates 
the major functional blocks of a direc- 
tion-insensitive analog tachometer that, 
except for its lack of directionality, offers 
the same advantages as does the analog 
tachometer described above. However, 
this circuit is conceptually simpler in 
that it does not contain integrators. 

This circuit processes the back-emf 
outputs of the two motor phases into a 
voltage directly proportional to magni- 
tude of the instantaneous rate of rota- 
tion of the shaft. As in the circuit de- 
scribed above, the processing in this 
circuit effects a straightforward combi- 
nation of mathematical operations 


leading to a final operation based on 
the identity (sin x) 2 + (cos x) 2 = 1 for 
any value of x. 

Further as in the circuit described 
above, the two back-emf signals are k 0sin 
6 for phase A and k 6 c os 6 for phase B, 
where k is a constant that depends on the 
electromagnetic characteristics of the 
motor. In the present case, the quantity 
that one seeks to measure is 1 6 1 . 

Each back-emf signal is fed to a dedi- 
cated squaring circuit. The outputs of 
the squaring circuits for phases A and B 
are thus proportional to (0sin 6) 2 and 
(0cos 6) 2 . The outputs of the squaring 
circuits are fed to an adder. By virtue of 
the identity (sin 6) 2 + (cos 6) 2 = 1 the 
output of the adder is proportional to 
6 2 ; this output is fed to a square-root cir- 
cuit to obtain a final output propor- 
tional to I 6 I . 

This work was done by David E. Howard 
and Dennis A. Smith of Marshall Space 
Flight Center. Further information is con- 
tained in a TSP ( seepage 1 ). 

This invention has been patented by NASA 
(U.S. Patent No. 6,084,398). Inquiries con- 
cerning nonexclusive or exclusive license for its 
commercial development should be addressed to 
Sammy Nabors, MSFC Commercialization As- 
sistance Lead, at sammy.a.nabors@nasa.gov. 
Refer to MFS-3 1142/3. 
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Software 


@ Algorithm-Based Fault Tol- 
erance for Numerical Sub- 
routines 

A software library implements a new 
methodology of detecting faults in nu- 
merical subroutines, thus enabling appli- 
cation programs that contain the subrou- 
tines to recover transparently from 
single-event upsets. The software library 
in question is fault-detecting middleware 
that is wrapped around the numerical- 
subroutines. Conventional serial versions 
(based on LAPACK and FFTW) and a 
parallel version (based on ScaLAPACK) 
exist. The source code of the application 
program that contains the numerical 
subroutines is not modified, and the 
middleware is transparent to the user. 

The methodology used is a type of algo- 
rithm-based fault tolerance (ABFT). In 
ABFT, a checksum is computed before a 
computation and compared with the 
checksum of the computational result; an 
error is declared if the difference between 
the checksums exceeds some threshold. 
Novel normalization methods are used in 
the checksum comparison to ensure cor- 
rect fault detections independent of algo- 
rithm inputs. In tests of this software re- 
ported in the peer-reviewed literature, this 
library was shown to enable detection of 
99.9 percent of significant faults while gen- 
erating no false alarms. 

This program was written by Michael Tur- 
mon, Robert Granat, and John Lou of Cal- 
tech for NASA’s Jet Propulsion Laboratory. 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-42193. 


Computational Support for 
Technology-Investment De- 
cisions 

Strategic Assessment of Risk and Tech- 
nology (START) is a user-friendly com- 
puter program that assists human man- 
agers in making decisions regarding 
research-and-development investment 
portfolios in the presence of uncertainties 
and of non-technological constraints that 
include budgetary and time limits, restric- 
tions related to infrastructure, and pro- 
grammatic and institutional priorities. 
START facilitates quantitative analysis of 
technologies, capabilities, missions, sce- 
narios and programs, and thereby enables 


the selection and scheduling of value-op- 
timal development efforts. START incor- 
porates features that, variously, perform 
or support a unique combination of func- 
tions, most of which are not systematically 
performed or supported by prior deci- 
sion-support software. These functions in- 
clude the following: 

• Optimal portfolio selection using an 
expected-utility-based assessment of ca- 
pabilities and technologies; 

• Temporal investment recommendations; 

• Distinctions between enhancing and 
enabling capabilities; 

• Analysis of partial funding for enhanc- 
ing capabilities; and 

• Sensitivity and uncertainty analysis. 
START can run on almost any comput- 
ing hardware, within Linux and related 
operating systems that include Mac OS X 
versions 10.3 and later, and can run in 
Windows under the Cygwin environment. 
START can be distributed in binary code 
form. START calls, as external libraries, 
several open-source software packages. 
Output is in Excel (.xls) file format. 

This program was written by Virgil Adu- 
mitroaie, Hook Hua, William Lincoln, Gary 
Block, Joseph Mrozinski, Kacie Shelton, 
Charles Weisbin, Alberto Elfes, and Jeffrey 
Smith of Caltech for NASA’s Jet Propulsion 
Laboratory. Further information is contained 
in a TSP (see page 1 ). 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-44800. 


DSN Resource Scheduling 

TIGRAS is client-side software, which 
provides tracking-station equipment 
planning, allocation, and scheduling 
services to the DSMS (Deep Space Mis- 
sion System) . TIGRAS provides func- 
tions for schedulers to coordinate the 
DSN (Deep Space Network) antenna 
usage time and to resolve the resource 
usage conflicts among tracking passes, 
antenna calibrations, maintenance, and 
system testing activities. TIGRAS pro- 
vides a fully integrated multi-pane 
graphical user interface for all schedul- 
ing operations. This is a great improve- 
ment over the legacy VAX VMS com- 
mand line user interface. 

TIGRAS has the capability to handle all 
DSN resource scheduling aspects from 
long-range to real time. TIGRAS assists 


NASA mission operations for DSN track- 
ing of station equipment resource re- 
quest processes from long-range load 
forecasts (ten years or longer), to mid- 
range, short-range, and real-time (less 
than one week) emergency tracking plan 
changes. TIGRAS can be operated by 
NASA mission operations worldwide to 
make schedule requests for the DSN sta- 
tion equipment. 

This program was written by Yeou-Fang 
Wang and John Baldwin of Caltech for 
NASAs Jet Propulsion Laboratory. 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-43771. 


tjl Distributed Operations 
Planning 

Maestro software provides a secure and 
distributed mission planning system for 
long-term missions in general, and the 
Mars Exploration Rover Mission (MER) 
specifically. Maestro, the successor to the 
Science Activity Planner, has a heavy em- 
phasis on portability and distributed op- 
erations, and requires no data replication 
or expensive hardware, instead relying on 
a set of services functioning on JPL insti- 
tutional servers. 

Maestro works on most current com- 
puters with network connections, includ- 
ing laptops. When browsing down-link 
data from a spacecraft, Maestro func- 
tions similarly to being on a Web 
browser. After authenticating the user, it 
connects to a database server to query an 
index of data products. It then contacts a 
Web server to download and display the 
actual data products. The software also 
includes collaboration support based 
upon a highly reliable messaging system. 
Modifications made to targets in one in- 
stance are quickly and securely transmit- 
ted to other instances of Maestro. 

The back end that has been developed 
for Maestro could benefit many future 
missions by reducing the cost of central- 
ized operations system architecture. 

This program was written by Jason Fox, 
Jeffrey Norris, Mark Powell, Kenneth Robe, 
and Khawaja Shams of Caltech for NASA’s 
Jet Propulsion Laboratory. 

This software is available for commercial 
licensing. Please contact Karina Edmonds of 
the California Institute of Technology at 
(626) 395-2322. Refer to NPO-43278. 
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Mechanics/Machinery 



# Phase-Oriented Gear Systems 

Larger mechanical advantages can be realized in smaller packages. 

Goddard Space Flight Center, Greenbelt, Maryland 




-Output Internal 
Ring Gear 


-Ground Internal 
Ring Gear 


Different Numbers of Teeth on the upper and lower part of each planet gear are used to accommo- 
date the phase differences occasioned by the use of different numbers of teeth on the ground and 
output gears. 


Phase-oriented gear systems are dif- 
ferential planetary transmissions in 
which each planet gear has two sets of 
unequal numbers of teeth indexed at 
prescribed relative angles (phases). 
The figure illustrates an application of 
the phase-oriented gearing concept to a 
relatively simple speed-reducing differ- 
ential planetary transmission that in- 
cludes a sun gear, an idler gear, three 
identical planet gears, a ground inter- 
nal ring gear, and an output internal 
ring gear. Typically, the ground internal 
ring gear and output internal ring gear 
have different numbers of teeth, giving 
rise to a progressive and periodic phase 
shift between the corresponding pairs 
of teeth engaged by each successive 
planet gear. To accommodate this 
phase shift, it is necessary to introduce 
a compensating phase shift between the 
ground-gear-engaging and output-gear- 
engaging sections of each planet gear. 
This is done by individually orienting 
each planet gear. 

Each planet gear consists of a top (out- 
put-gear-engaging) part and a bottom 
(ground-gear-engaging) part that are 
coaxially joined into a unit. Typically, the 
numbers of teeth in the top and bottom 
parts differ by one, (but could differ by 
more than one). Typically, the teeth on 
the bottom part are indexed with those 
on the top part such that at least one 
tooth on the bottom part lies at the same 
angular position as does one tooth on the 
top part; in other words, these two teeth 
are in phase with each other. The next 
tooth on the top part is slightly out of 
phase with the nearest corresponding 
tooth on the bottom part. Proceeding cir- 
cumferentially around the planet gear, 
this phase difference increases linearly, 
tooth-by-tooth, until the total phase 
change reaches 360° upon completing 
one full rotation. Depending on the spe- 
cific numbers of teeth, such an arrange- 
ment could offer many tooth-phasing 
possibilities, making it possible to include 
various numbers of appropriately phased 
(oriented) planet gears. This is advanta- 
geous because in some applications, it is 
desirable to increase numbers of planet 


gears in order to make more teeth share 
the total load, thereby reducing loads on 
individual teeth. 

The phase-oriented gearing concept 
admits of so many variations that is not 
possible to describe or even to merely 


list them within the space available for 
this article. In general, it can be said 
that the numbers of teeth and the an- 
gles can be chosen to achieve desired 
speed ratios and other design goals. The 
phase-oriented gearing concept can be 
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implemented in conjunction with other 
advanced planetary-gearing concepts, 
including those of carrierless and gear- 
bearing transmission types described in 
several NASA Tech Briefs articles in re- 
cent years. Relative to prior differential 
planetary transmissions, phase-oriented 
gear systems offer advantages of simplic- 
ity, ruggedness, strength, and smooth- 


ness of operation. Moreover, for a typi- 
cal speed-reducing transmission, phase- 
oriented gearing makes it possible to 
obtain unprecedentedly large mechani- 
cal advantage in a compact package — 
in effect, what previously would have 
been characterized as two-stage per- 
formance in the space previously occu- 
pied by a one-stage package. 


This work was done by John M. Vranish 
of Goddard Space Flight Center. 

This invention is owned by NASA, and a 
patent application has been filed. Inquiries 
concerning nonexclusive or exclusive license 
for its commercial development should be ad- 
dressed to the Patent Counsel, Goddard Space 
Flight Center, (301 ) 286-7351. Refer to GSC- 
14790-1. 
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Manufacturing & Prototyping 


i=! Freeze Tape Casting of Functionally Graded Porous Ceramics 

Pore structures can be tailored in ways heretofore impossible. 

John H. Glenn Research Center, Cleveland, Ohio 


Freeze tape casting is a means of mak- 
ing preforms of ceramic sheets that, 
upon subsequent completion of fabrica- 
tion processing, can have anisotropic 
and/or functionally graded properties 
that notably include aligned and graded 
porosity. Freeze tape casting was devel- 
oped to enable optimization of the mi- 
crostructures of porous ceramic compo- 
nents for use as solid oxide electrodes in 
fuel cells: Through alignment and grad- 
ing of pores, one can tailor surface areas 
and diffusion channels for flows of gas 
and liquid species involved in fuel-cell 
reactions. Freeze tape casting offers sim- 
ilar benefits for fabrication of optimally 
porous ceramics for use as catalysts, gas 
sensors, and filters. 

Freeze tape casting includes, and goes 
beyond, traditional tape casting, in 
which an aqueous ceramic slip (ceramic 
and/or ceramic precursor particles sus- 
pended in water) is cast onto a polyeth- 
ylene terephthalate) or poly(tetrafluo- 
roethylene) carrier film by use of a 
doctor-blade assembly. The slip may also 


contain one or more organic solvent (s) 
plus a significant quantity of organic 
binders that make the tape strong and 
flexible after evaporation of the sol- 
vents). Traditional tape casting has 
been used to make ceramic sheets, for 
electronic and structural applications, 
ranging in thickness from 5 pm to 1,000 
pm. Freeze tape casting expands this 
range to greater than 3 mm, thus elimi- 
nating the need for lamination steps in 
traditional processing. 

An apparatus for freeze tape casting 
(see figure) is basically a traditional 
tape-casting apparatus augmented with a 
freezing bed. By use of the freezing bed, 
the water in the slip is frozen, from the 
bottom up through the thickness, imme- 
diately after the tape has been cast. For 
fabricating a porous ceramic, the freez- 
ing affords several advantages over tradi- 
tional tape casting: 

• In traditional tape casting, the water and 
any organic solvents are allowed to evap- 
orate slowly before the tape is sintered 
to form the final ceramic. The evapora- 


tion can cause or be accompanied by 
undesired compositional and physical 
changes, including settling of particles 
out of suspension and density gradients. 
The freezing prevents such changes. 

• Water alone or high melting point or- 
ganic solvents (including benzene, cy- 
clohexane, tertiary butyl alcohol, and 
camphene), can become formed into 
channels or other unique pore struc- 
tures when directionally solidified. 
The channels or other pore structures 
can be tailored to some extent 
through control of the rate of freezing 
and/or the addition of such freezing 
additives as glycerol, glycols, anti- 
freeze proteins, and/or alcohols. 

• Porosity can also be tailored through 
choice of the concentration of sus- 
pended ceramic particles: typically, 
freeze-tape-cast ceramics contain open 
pore structures when the slips contain 
less than 45 volume percent of ceramic 
and/or metallic solids. 

• In freeze tape casting, the solids loading 
can be made as low as 5 volume percent 


Doctor-Blade 

Assembly 



FREEZE-TAPE-CASTING APPARATUS 


Room 
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In Freeze Tape Casting, a slip is cast into a tape as in traditional tape casting, but then, unlike traditional tape casting, the water and any other solvents in the tape 
are frozen. The frozen tape is then freeze-dried before sintering. (Temperature profile is indicated for an aqueous system) 
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to obtain exceptionally high porosity. 

After casting and freezing, the frozen 
tape is diced into sections for freeze- 
drying, in which the water and any 
other solvents are removed by sublima- 
tion. Because, in sublimation, the solid- 
ified solvents are transformed into 
gases without passing through interme- 
diate liquid phases, there are no capil- 
lary forces like those associated with liq- 
uid-to-vapor transitions that occur 


during drying in traditional tape cast- 
ing. Because of the absence of capillary 
forces, the changes in the dimensions 
of microstructures and in the overall 
thickness of the tape are negligible and 
hence processing is simplified with the 
absence of cracks and other drying de- 
fects notable in traditional tapes. After 
freeze-drying, the sections of tape are 
cut into the desired shape, then sin- 
tered. 


This work was done by Stephen W. Sofie of 
QSS Group, Inc. for Glenn Research Center. 
Further information is contained in a TSP 
( see page 1 ). 

Inquiries concerning rights for the com- 
mercial use of this invention should be ad- 
dressed to NASA Glenn Research Center, In- 
novative Partnerships Office, Attn: Steve 
Fedor, Mail Stop 4-8, 21000 Brookpark 
Road, Cleveland, Ohio 44135. Refer to 
LEW-1 7 628-1. 


§ Electrophoretic Deposition on Porous Non-Conductors 

EPD is simplified and made more widely applicable. 

John H. Glenn Research Center, Cleveland, Ohio 



Particles of YSZ Were Deposited on an NiO-YSZ substrate in this experimental EPD apparatus. 


A method of electrophoretic deposi- 
tion (EPD) on substrates that are porous 
and electrically non-conductive has 
been invented. Heretofore, in order to 
perform an EPD, it has been necessary 
to either ( 1 ) use a substrate material that 
is inherently electrically conductive or 
(2) subject a non-conductive substrate to 
a thermal and/ or chemical treatment to 
render it conductive. 

In the present method, instead of rely- 
ing on the electrical conductivity of the 
substrate, one ensures that the substrate 
is porous enough that when it is im- 
mersed in an EPD bath, the solvent pen- 
etrates throughout the thickness, thereby 
forming quasi-conductive paths through 
the substrate. By making it unnecessary 
to use a conductive substrate, this 
method simplifies the overall EPD 
process and makes new applications pos- 
sible. The method is expected to be espe- 
cially beneficial in enabling deposition of 
layers of ceramic and/ or metal for chem- 
ical and electrochemical devices, notably 
including solid oxide fuel cells. 

In the initial application for which this 
method was invented, there is a require- 
ment to fabricate unitary anode/solid 
electrolyte structures for solid oxide fuel 
cells, each such structure comprising the 
following: 

• A porous cermet anode layer made of 
nickel and yttria-stabilized zirconia 
(YSZ) and 

• A solid electrolyte in the form a thin 
dense layer of YSZ. 

The starting material for fabricating 
the structure is a porous layer consist- 
ing mostly of pellets of a composite of 
NiO and YSZ. The dense layer of YSZ is 
deposited on the NiO-YSZ composite 
later by means of EPD according to the 
present method, then the resulting 


(NiO-YSZ) /YSZ structure is sintered. 
The final fabrication step takes place 
after the sintered structure has been 
placed in a fuel cell: The reducing at- 
mosphere maintained on the anode 
side of the fuel cell during testing re- 
duces the NiO to metallic Ni, thereby 
creating the desired Ni anode having 
the porosity needed for transport of re- 
actant and product gases. 

The figure depicts the apparatus used 
in experiments to develop and demon- 
strate the present EPD method. The ap- 
paratus included a reservoir containing 
acetylacetone (as a nominally nonaque- 
ous solvent) . YSZ powder was stably sus- 
pended in the solvent by electrostatic re- 
pulsion, and YSZ particles were 
replenished to the depositing electrode 
by use of a magnetic stirrer. The YSZ par- 
ticles acquired positive surface electric 


charges in a manner similar to that of 
surface charging in an aqueous 
medium. This charging has been conjec- 
tured to be a consequence of the pres- 
ence of residual water in the acetylace- 
tone. Regardless of the exact physical 
mechanism, the charging was fortuitous 
in that it contributed to the desired elec- 
trophoretic effect. 

The NiO-YSZ deposition substrate 
was mounted on a stationary poly(tetra- 
fluoroethylene) electrode holder and 
connected via a spring contact to the 
negative side of a DC power supply so 
that it became a deposition electrode 
(cathode). A platinum counter elec- 
trode was mounted on a movable 
poly(tetrafluoroethylene) electrode 
holder and connected to the positive 
side of the power supply, so that it be- 
came an anode. For most experimental 
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depositions, the distance between the 
electrodes was kept at 10 mm. The dep- 
ositions were carried out at various con- 
stant applied voltages ranging from 50 
to 300 V for periods of 30 s to 5 min. 
During each deposition, the force of 
the applied electric field on suspended 
YSZ particles caused these particles to 
migrate towards the cathode and even- 


tually deposit. A solid-oxide fuel cell 
containing an anode/ solid-electrolyte 
structure fabricated in part by this 
method exhibited an areal power den- 
sity of 611 W/cm 2 in operation at a tem- 
perature of 800 °C. 

This work was done by Charles Compson, 
Laxmidhar Besra, and Meilin Liu of Georgia 
Institute of Technology for Glenn Research 


Center. Further information is contained in a 
TSP (seepage 1 ). 

Inquiries concerning rights for the com- 
mercial use of this invention should be ad- 
dressed to NASA Glenn Research Center, In- 
novative Partnerships Office, Attn: Steve 
Fedor, Mail Stop 4-8, 21000 Brookpark 
Road, Cleveland, Ohio 44135. Refer to 
LEW-18009-1. 
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Bio-Medical 


@Two Devices for Removing Sludge From Bioreactor 
Wastewater 

These devices can operate continuously and are self-cleaning. 

Lyndon B. Johnson Space Center, Houston, Texas 


Two devices — a magnetic separator and 
a special filter denoted a self-regenerating 
separator (SRS) — have been developed 
for separating sludge from the stream of 
wastewater from a bioreactor. These de- 
vices were originally intended for use in mi- 
crogravity, but have also been demon- 
strated to function in normal Earth gravity. 

The magnetic separator (see Figure 
1) includes a thin-walled nonmagnetic, 
stainless-steel cylindrical drum that ro- 
tates within a cylindrical housing. The 
wastewater enters the separator 


through a recirculation inlet, and 
about 80 percent of the wastewater flow 
leaves through a recirculation outlet. 
Inside the drum, a magnet holder posi- 
tions strong permanent magnets sta- 
tionary and, except near a recirculation 
outlet, close to the inner drum surface. 
To enable magnetic separation, mag- 
netite (a ferromagnetic and magneti- 
cally soft iron oxide) powder is mixed 
into the bioreactor wastewater. The 
magnetite becomes incorporated into 
the sludge by condensation, onto the 


powder particles, of microbe flocks that 
constitute the sludge. As a result, the 
magnets inside the drum magnetically 
attract the sludge onto the outer sur- 
face of the drum. 

As the drum rotates, it carries the at- 
tracted sludge to the recirculation out- 
let; the decrease in the magnetic field 
at this location releases the sludge 
from the drum, making it possible for 
the sludge to be flushed out by the re- 
circulation-outlet flow. A poly(tetraflu- 
oroethylene) scraper aids in removing 
any remnants of sludge that may still 
stick to the drum. An effluent outlet is 
positioned as far as possible from the 
recirculation ports in order to provide 
the longest possible separation path. 
The effluent flow amounts to about 20 
percent of the inlet flow. The flows of 
sludge and effluent are controlled by 
use of a recirculation pump and an ef- 
fluent pump, respectively. 

The SRS (see Figure 2) includes two 
concentric tubes with relatively coarse 
filtering material (pipe cleaners) in the 
annular space between them. The 
wastewater is recirculated through the 
inner tube. Through an axial slot in 
the inner tube, a portion of the flow is 
diverted into the annular space be- 
tween the tubes. This portion of the 
flow becomes filtered in the annular 
space and is drawn off as effluent 
through an axial slot in the outer tube 
on the side opposite that of the slot in 
the inner tube. 

Eventually, any conventional filter be- 
comes clogged and thereby loses hy- 
draulic capacity and must be replaced. 
An ideal filter would have to continuously 
clean itself so as not to impede the flow of 
the filtrate. The SRS acts as such an ideal 
or nearly ideal filter because of the self-re- 
generating aspect of its operation: 

• Initially, sludge builds up quickly in 
the pipe-cleaner flow field, where the 
sludge thus trapped serves as a fine fil- 
ter medium for incoming sludge parti- 
cles. As more sludge builds up toward 
the interior of the inner tube, it is ex- 


Recirculation Outlet 


t 



Figure 1. The Magnetic Separator rapidly separates magnetic sludge from wastewater. 
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posed to the fast recirculation flow, 
which washes away newly settled 
sludge. A dynamic equilibrium is 
quickly reached between the deposi- 
tion and removal of sludge. 

• Unlike what would happen in a conven- 
tional filter, the buildup of sludge in the 
pipe-cleaner flow field does not con- 
tinue indefinitely until excessive clog- 
ging occurs. Instead, the SRS design uti- 
lizes the fact that sludge is a living 
organism that digests its own dead parts 
stagnating in the pipe-cleaner flow 


field. The digestion process continu- 
ously opens up new filter pores that 
make it possible for water to continue 
to flow. This process gives rise to a sec- 
ond dynamic equilibrium — between 
clogging and the generation of new 
flow paths in the pipe-cleaner flow field. 

Only the first pipe cleaner is needed as 
a filter element; the other pipe cleaners 
are needed only for support of the tubes 
in their coaxial alignment. Prior to test- 
ing, it was expected that sludge would 
slowly migrate through the pipe-cleaner 


flow field and reach the effluent port, 
where it would reduce the quality of the 
filtrate. Instead, it was found that the 
sludge remained stationary at the first 
pipe cleaner and no migration was ob- 
served, even after two weeks of operation. 

This work was done by Shivaun Archer, G. 
Duncan Hitchens, Harry Jabs, Jennifer 
Cross, Michelle Pilkinton, and Michael Tay- 
lor of Lynntech, Inc. for Johnson Space Cen- 
ter. For further information, contact the John- 
son Innovative Partnerships Office at (281) 
483-3809. Refer to MSC-23293 


@ Portable Unit for Metabolic Analysis 

Respiratory signals can be temporally resolved within respiratory cycles. 

John H. Glenn Research Center, Cleveland, Ohio 


The Portable Unit for Metabolic Analy- 
sis (PUMA) is an instrument that meas- 
ures several quantities indicative of 
human metabolic function. Specifically, 
this instrument makes time-resolved 
measurements of temperature, pressure, 
flow, and the partial pressures of oxygen 
and carbon dioxide in breath during both 
inhalation and exhalation. 

Portable instruments for measuring 
these quantities have been commercially 
available, but the response times of those 
instruments are too long to enable tem- 
poral resolution of phenomena on the 
time scales of human respiration cycles. 
In contrast, the response time of the 
PUMA is significantly shorter than char- 
acteristic times of human respiration phe- 
nomena, making it possible to analyze 
varying metabolic parameters, not only 
on sequential breath cycles but also at 
successive phases of inhalation and exha- 
lation within the same breath cycle. 

In operation, the PUMA is positioned 
to sample breath near the subject’s 
mouth. Commercial off-the-shelf sensors 


are used for three of the measurements: a 
miniature pressure transducer for pres- 
sure, a thermistor for temperature, and 
an ultrasonic sensor for flow. Sensors de- 
veloped at Glenn Research Center are 
used for measuring the partial pressures 
of oxygen and carbon dioxide: 

• The carbon dioxide sensor exploits 
the relatively strong absorption of in- 
frared light by carbon dioxide. Light 
from an infrared source passes 
through the stream of inhaled or ex- 
haled gas and is focused on an in- 
frared-sensitive photodetector. 

• The oxygen sensor exploits the effect 
of oxygen in quenching the fluores- 
cence of ruthenium-doped organic 
molecules in a dye on the tip of an op- 
tical fiber. A blue laser diode is used to 
excite the fluorescence, and the opti- 
cal fiber carries the fluorescent light to 
a photodiode, the temporal variation 
of the output of which bears a known 
relationship with the rate of quench- 
ing of fluorescence and, hence, with 
the partial pressure of oxygen. 


The outputs of the sensors are digi- 
tized, preprocessed by a small onboard 
computer, and then sent wirelessly to a 
desktop computer, where the collected 
data are analyzed and displayed. In addi- 
tion to the raw data on temperature, 
pressure, flow, and mole fractions of oxy- 
gen and carbon dioxide, the display can 
include volumetric oxygen consump- 
tion, volumetric carbon dioxide produc- 
tion, respiratory equivalent ratio, and 
volumetric flow rate of exhaled gas. 

This work was done by Daniel L. Dietrich, 
Nancy D. Piltch, Mark E. Lewis, Jeffrey R. 
Juergens, Michael J. Lichter, Peter M. Struk, 
and Dale M. Diedrick of Glenn Research 
Center; Russell W. Valentine of Case West- 
ern Reserve University; and Richard D. Pet- 
tegrew of the National Center for Micrograv- 
ity Research. 

Inquiries concerning rights for the commer- 
cial use of this invention should be addressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21 000 Brookpark Road, Cleveland, 
Ohio 44135. Refer to LEW-17945-1. 
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# Flash Diffusivity Technique Applied to Individual Fibers 

Inversion of spatially and temporally resolved temperature measurements yields 
thermal-diffusivity maps. 

Lyndon B. Johnson Space Center, Houston, Texas 


A variant of the flash diffusivity tech- 
nique has been devised for determining 
the thermal diffusivities, and thus the 
thermal conductivities, of individual 
aligned fibers. The technique is intended 
especially for application to nanocompos- 
ite fibers, made from narrower fibers of 
polyphenylene benzobisthiazole (PBZT) 
and carbon nanotubes. These highly 
aligned nanocomposite fibers could ex- 
ploit the high thermal conductivities of 
carbon nanotubes for thermal-manage- 
ment applications. 

In the flash diffusivity technique as 
practiced heretofore, one or more heat 
pulse (s) is (are) applied to the front face 
of a plate or disk material specimen and 
the resulting time-varying temperature 
on the rear face is measured. Usually, the 


heat pulse is generated by use of a xenon 
flash lamp, and the variation of tempera- 
ture on the rear face is measured by use 
of an infrared detector. The flash energy 
is made large enough to produce a use- 
fully high temperature rise on the rear 
face, but not so large as to significantly 
alter the specimen material. Once the 
measurement has been completed, the 
thermal diffusivity of the specimen is 
computed from the thickness of the spec- 
imen and the time dependence of the 
temperature variation on the rear face. 

Heretofore, the infrared detector 
used in the flash diffusivity technique 
has been a single-point detector, which 
responds to a spatial average of the ther- 
mal radiation from the rear specimen 
surface. Such a detector cannot distin- 


guish among regions of differing diffu- 
sivity within the specimen. Moreover, 
two basic assumptions of the thermal- 
diffusivity technique as practiced 
heretofore are that the specimen is ho- 
mogeneous and that heat flows one-di- 
mensionally from the front to the rear 
face. These assumptions are not valid 
for an inhomogeneous (composite) ma- 
terial. 

In the present variant of the flash dif- 
fusivity technique, one uses an infrared 
electronic camera fitted with a micro- 
scope lens to record the spatial as well as 
the temporal variations in thermal radia- 
tion emitted from the rear face of the 
specimen (see figure). In the recorded 
image of a composite-material speci- 
men, it is possible to distinguish between 



An Infrared Camera With Microscope Lens measures the spatial as well as temporal variations of temperature on the rear face of a specimen that has been 
flash-heated on its front face. 
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individual fibers, or between a fiber and 
the surrounding matrix material. 
Hence, it is possible to measure the rear- 
face temperature variations of individual 
fibers. These variations can be corre- 
lated with predictions of a computa- 
tional model of heat transfer in the com- 
posite specimen to obtain a diffusivity 
map of the specimen. 

The technique was demonstrated on 
a specimen containing pure PBZT 
fibers, one nanocomposite PBZT/car- 
bon-nanotube fiber, and one copper 


fiber mounted longitudinally in an 
epoxy matrix. The copper fiber, having 
known thermal conductivity, was in- 
cluded for qualitative comparison. The 
temperature transients of the pure 
PBZT fibers were not distinguishable 
from that of the matrix, and the ther- 
mal diffusivity of the matrix and PBZT 
fibers was found to be 0.0032 cm 2 /s. 
The thermal diffusivity of the 
PBZT/carbon-nanotube composite 
fiber was found to be 0.049 cm 2 /s; the 
true bulk diffusivity of the PBZT/car- 


bon-nanotube composite could be 
higher than the value computed from 
the measurements because the heat- 
transfer model used in the computa- 
tions does not account for thermal cou- 
pling between the fibers and the matrix. 

This work was done by Brian Mayeaux, 
Leonard Yowell of Johnson Space Center, and 
Hsin Wang of Oak Ridge National Labora- 
tory. For further information, contact the 
Johnson Innovative Partnerships Office at 
(281) 483-3809. 

MSC-23609 


9 System for Thermal Imaging of Hot Moving Objects 

This advanced ratio-intensity thermography system has numerous potential uses. 

Langley Research Center, Hampton, Virginia 


Filter 

Wheel 


Barlow 

Lens 




Near-Infrared 

Camera 



Infrared Radiation 
From Observed 
Object 


Telescope 


The RITS Acquires Images at Multiple Wavelengths in rapid succession by means of synchronized op- 
eration of a filter wheel and near-infrared camera operating at a frame rate of 30 Hz. The design and 
mode of operation of this system are such that if two wavelengths are used, then a total of 15 pairs 
of complete images is acquired during 1 second. 


The High Altitude/Re-Entry Vehicle 
Infrared Imaging (HARVTI) system is a 
portable instrumentation system for 
tracking and thermal imaging of a possi- 
bly distant and moving object. The 
HARVII is designed specifically for 
measuring the changing temperature 
distribution on a space shuttle as it reen- 
ters the atmosphere. The HARVTI sys- 
tem or other systems based on the de- 
sign of the HARVTI system could also be 
used for such purposes as determining 
temperature distributions in fires, on 
volcanoes, and on surfaces of hot mod- 
els in wind tunnels. In yet another po- 
tential application, the HARVTI or a sim- 
ilar system would be used to infer 
atmospheric pollution levels from im- 
ages of the Sun acquired at multiple 
wavelengths over regions of interest. 

The HARVII system includes the 
Ratio Intensity Thermography System 
(RITS) and a tracking subsystem that 
keeps the RITS aimed at the moving ob- 
ject of interest. The subsystem of pri- 
mary interest here is the RITS (see fig- 
ure), which acquires and digitizes 
images of the same scene at different 
wavelengths in rapid succession. Assum- 
ing that the time interval between suc- 
cessive measurements is short enough 
that temperatures do not change appre- 
ciably, the digitized image data at the 
different wavelengths are processed to 
extract temperatures according to the 
principle of ratio-intensity thermogra- 
phy: The temperature at a given location 
in a scene is inferred from the ratios be- 
tween or among intensities of infrared 
radiation from that location at two or 
more wavelengths. This principle, based 
on the Stefan-Boltzmann equation for 


the intensity of electromagnetic radia- 
tion as a function of wavelength and 
temperature, is valid as long as the ob- 
served body is a gray or black body and 
there is minimal atmospheric absorp- 
tion of radiation. 

In this particular variant of ratio-in- 
tensity thermography, the wavelengths 
are chosen to be in the near infrared 
part of the spectrum and are selected to 
minimize effects of atmospheric absorp- 
tion. A nearly ideal pair of wavelengths 
for minimizing effects of atmospheric 
absorption would be 1.3 and 1.69 pm. 
An additional advantage of this choice 
of wavelengths is that conventional, 
high-quality, commercial optics work 
well at wavelengths below about 2 pm; if 
longer wavelengths were chosen, it 
would be necessary to design and build 
special infrared optical components. 

In practice, it was necessary to choose 
1.55 pm as the longer of the two wave- 
lengths because the infrared camera 


available at the time of development of 
the RITS was found not to be usable out 
to 1.69 pm. In addition, the pass bands 
of the filters used to select these observa- 
tion wavelengths were required to be 10 
to 20 nm wide in order to enable the ac- 
quisition of sufficiently strong signals. 
The finite widths of these pass bands 
contribute to variations of atmospheric 
effects at different water-vapor levels. 
The use of three or more wavelengths 
could enable almost complete correc- 
tion for these effects. Alternately, if one 
continued to use two wavelengths, then 
corrections could be derived from meas- 
urements made by a spectrometer 
aimed at the Sun. 

This work was done by Leonard Weinstein 
and Jason Hundley of Langley Research Cen- 
ter. For more information, download the 
Technical Support Package (free white 
paper) at www.techhriefs.com/tsp under 
the Physical Sciences category. 

LAR- 16385-1 
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Large Solar-Rejection Filter 

This lightweight filter comprises a multilayer spectral coating on a flexible 
membrane. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An optical filter consisting of a multi- 
layer spectral coating on a flexible mem- 
brane has been designed to be placed in 
front of the 200-in. (5.08-m) Hale tele- 
scope on Mt. Palomar. The filter is in- 
tended to protect the telescope against 
solar radiant flux and limit solar heating 
of the interior of the telescope dome 
while transmitting light at the 1,064-nm 
wavelength of the Mars Laser Communi- 
cation Demonstration. 

For supporting a multilayer spectral 
coating in this application, the flexible 
membrane was chosen as a lightweight, 
less-expensive alternative to a conven- 
tional thick optical-glass substrate. Multi- 
layer spectral coatings have been used 
for decades, and membranes have more 
recently come into use as substrates for 
mirrors. However, until now, there have 
been few (if any) published instances of 


multilayer coating of membranes to 
form lightweight solar-rejection/narrow- 
band-pass filters. 

The main problem to be solved in de- 
signing the present filter was to satisfy 
both (1) the need for a large number 
(hundreds in a typical first approxima- 
tion) of coating layers needed to obtain 
the desired broad-band-solar-rejection/ 
narrow-band-pass spectral characteristic 
and (2) the need to limit the number of 
layers to no more than the maximum 
(~60) that the membrane could support. 
The solution is a design of fewer than 50 
layers having the following features: 

• The front surface of the membrane is 
coated with 25 dielectric layers alter- 
nating between higher and lower in- 
dices of refraction. These layers are de- 
signed to efficiently reflect the visible 
and near-infrared (wavelengths up to 


1.0 pm) light in which the Sun pre- 
dominantly radiates. 

• The back surface of the membrane is 
coated with very thin layers of copper 
and silicon chosen to reject the re- 
mainder of the spectrum. 

• Multiple dielectric layers stacks are de- 
posited on what would otherwise be 
the exposed opposite outer surfaces of 
the aforementioned layers. The mate- 
rials and thicknesses of these layers are 
chosen to induce the desired narrow 
passband, centered at 1,064 nm, in the 
broad-band-reject spectrum of the sili- 
con and copper layers. 

This work was done by William Roberts of 
Caltech, David Sheikh of Surface Optics Corpo- 
ration, and Brian Patrick of SRS Corporation 
for NASA’s Jet Propulsion Laboratory. For more 
information, contact iaoffice@jplnasa.gov 
NPO-41942 


® Improved Readout Scheme for SQUID-Based Thermometry 

Flux-unit-counting ambiguities would be eliminated. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An improved readout scheme has 
been proposed for high-resolution ther- 
mometers, (HRTs) based on the use of 
superconducting quantum interference 
devices (SQUIDs) to measure tempera- 
ture-dependent magnetic susceptibili- 
ties. The proposed scheme would elimi- 
nate counting ambiguities that arise in 
the conventional scheme, while main- 
taining the superior magnetic-flux sensi- 
tivity of the conventional scheme. The 
proposed scheme is expected to be espe- 
cially beneficial for HRT-based tempera- 
ture control of multiplexed SQUID- 
based bolometer sensor arrays. 

SQUID-based HRTs have become 
standard for measuring and controlling 


temperatures in the sub-nano-Kelvin 
temperature range in a broad range of 
low-temperature scientific and engi- 
neering applications. A typical SQUID- 
based HRT that utilizes the conven- 
tional scheme includes a coil wound on 
a core made of a material that has tem- 
perature-dependent magnetic suscepti- 
bility in the temperature range of inter- 
est. The core and the coil are placed in 
a DC magnetic field provided either by 
a permanent magnet or as magnetic 
flux inside a superconducting outer 
wall. The aforementioned coil is con- 
nected to an input coil of a SQUID. 
Changes in temperature lead to 
changes in the susceptibility of the core 


and to changes in the magnetic flux de- 
tected by the SQUID. 

The SQUID readout instrumentation 
is capable of measuring magnetic-flux 
changes that correspond to temperature 
changes down to a noise limit ~0.1 
nK/Hz 1/2 . When the flux exceeds a few 
fundamental flux units, which typically 
corresponds to a temperature of -100 
nK, the SQUID is reset. The tempera- 
ture range can be greatly expanded if 
the reset events are carefully tracked and 
counted, either by a computer running 
appropriate software or by a dedicated 
piece of hardware. 

While adequate for many applica- 
tions, the conventional scheme has 
drawbacks: If the temperature is 

changed rapidly or the temperature 
noise is high, the counting hardware 
and/ or software loses flux count. In the 
case of a software counter, the tempera- 
ture reading is lost entirely if the soft- 
ware is reset or restarted. In the case of 
a multiplexed SQUID controller, these 
drawbacks become more severe because 
flux readings are taken less frequently. 



A Branch Containing a Secondary SQUID would be added to a single-SQUID HRT circuit. The second- 
ary SQUID would register a small fraction of the flux passing through the primary SQUID. 
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The proposed scheme is intended to 
eliminate these drawbacks. The scheme 
calls for including a secondary SQUID 
and its readout instrumentation that 
would register a small fraction of the 
magnetic flux passing through a primary 
SQUID. The scheme includes the fol- 
lowing elements: 

• Winding a secondary coil of fewer 
turns around the core to a second 
readout; or 

• Forming a circuit branch parallel to 
the main coil with the secondary 
SQUID input coil in series with a large 
(compared to the SQUID input coil in- 
ductance) inductor; or 


• Forming a circuit branch parallel to 
the main coil with a large inductor in 
series with a SQUID input coil shunted 
by a small inductor (see figure). 

The goal is to avoid having to reset 
the secondary SQUID in the tempera- 
ture range of interest, while maintain- 
ing the capability of determining the 
flux state of the primary SQUID unam- 
biguously. If the secondary SQUID 
readout were monitored by a 16-bit 
data-acquisition board and the digitiza- 
tion effected by the board determined 
the readout accuracy, then the dynamic 
range afforded by this scheme could be 
optimized by designing the secondary 


SQUID readout to be about 1/32,000 
as sensitive as is the primary SQUID 
readout. In a typical application, this 
level of secondary-SQUID sensitivity 
would correspond to a temperature 
range ~3 mK. In this temperature 
range, there would be no need to ac- 
tively track the flux state to maintain fi- 
delity of the readout. To avoid the need 
for counting hardware altogether, a ter- 
tiary readout could be added. 

This work was done by Konstantin Penanen 
of Caltech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page 1 ). 

NPO-41757 
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©Error Rates and Channel Capacities in Multipulse PPM 

It is now possible to compare expected performances of candidate modulation schemes. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


A method of computing channel ca- 
pacities and error rates in multipulse 
pulse-position modulation (multipulse 
PPM) has been developed. The method 
makes it possible, when designing an op- 
tical PPM communication system, to de- 
termine whether and under what condi- 
tions a given multipulse PPM scheme 
would be more or less advantageous, rel- 
ative to other candidate modulation 
schemes. 

In conventional M - ary PPM, each sym- 
bol is transmitted in a time frame that is 
divided into M time slots (where M is an 
integer >1) , defining an M-symbol alpha- 
bet. A symbol is represented by transmit- 
ting a pulse (representing “1”) during 
one of the time slots and no pulse (rep- 
resenting “0”) during the other M - 1 
time slots. Multipulse PPM is a general- 
ization of PPM in which pulses are trans- 
mitted during two or more of the M time 
slots. If the number of pulses per symbol 
is n , then the number of symbols in the 


alphabet is given by the binomial coeffi- 
cient 

C™ =M\/[n\(M-n)\]. 

The method is based partly on an 
analysis of the conditional probability, 
pi(y) or po(y), that the actual value, y, 
of the noisy signal detected in a re- 
ceiver during a given time slot repre- 
sents a transmitted 1 or a transmitted 0, 
respectively. For purposes of the analy- 
sis, the signal-propagation channel is 
assumed to be memoryless. The analy- 
sis includes consideration of L(y) = 
pi(y)/po(y)> defined as the likelihood 
ratio for receiving value y during the 
time slot. It is assumed that L(y) is fi- 
nite and, as is true for many channels 
of practical interest, that L(y) increases 
monotonically with y. This analysis 
leads to a maximum-likelihood (ML) 
decision rule to be used by the receiver. 
The rule turns out to be equivalent to 



SER Values were calculated for 16-ary PPM using several different values n and two different noise lev- 
els in a Poisson channel. 


the intuitive conjecture that the symbol 
most likely to have been transmitted, 
given the actual signal values detected 
during a symbol frame, is the one rep- 
resented by 1 in each of the n time slots 
for which the detected signal values 
were the largest. 

Next, an analysis of the probabilities 
that ML decisions are correct or incor- 
rect for various detected signal levels 
leads to the following equation for the 
symbol error rate (SER) — that is, the 
probability of incorrect symbol detection: 

oo M—n n 
5=0 v=0 t = 1 

[l-P 1 (i) B -*]C"-"/ Ib (*) 0 Po(*-l) M_B_0 , 


where I(t,v) =1/ C t t+V , P A (s) is the proba- 
bility that the digitized value of the de- 
tected signal (e.g., the number of de- 
tected photons) in a “1” slot is less than 
or equal to 5, and P 0 (s) is the probability 
that the digitized value of the detected 
signal in a “0” slot is less than or equal to 
5 . The figure presents an example of 
SER values calculated by use of these 
equations. 

Next, a comparative analysis of 
throughput achievable in conventional 
and multipulse PPM under bandwidth, 
average-power, and peak-power con- 
straints leads to the following equation 
for the channel capacity: 


C = -E Y \ x=1 log 


pm 

P(Y) 


M — n PoW) i • / i 

~ e y |x=o log bits/slot, 


M 


where p(Y) = ( n/M)p l (Y ) + [(M- n)/M\ 
po(Y) is denoted the probability mass 
function for a randomly chosen slot and 
E Y \x=i or Ey\ x =o is the expected value of 
signal level Y in a slot for which the 
transmitted signal value, X, was 1 or 0, 
respectively. 

This work was done by Jon Hamkins and 
Bruce Moision of Caltech for NASA’s Jet 
Propulsion Laboratory. For further informa- 
tion, contact iaoffice@jpl.nasa.gov. 
NPO-41536 
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©Two Mathematical Models of Nonlinear Vibrations 

Model parameters are fit to empirical vibration data. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


Two innovative mathematical models 
of nonlinear vibrations, and methods of 
applying them, have been conceived as 
byproducts of an effort to develop a 
Kalman filter for highly precise estima- 
tion of bending motions of a large truss 
structure deployed in outer space from a 
space-shuttle payload bay. These models 
are also applicable to modeling and 
analysis of vibrations in other engineer- 
ing disciplines, on Earth as well as in 
outer space. 

The first model is denoted the ampli- 
tude-dependent stiffness (ADS) model to 
emphasize the difference between it and 
the classical linear harmonic-oscillator 
model, in which stiffness is a constant. The 
ADS model is embodied in the equation 

x + £x + K(x, x)x = 0, 

where x is the instantaneous amplitude 
of the oscillating position or modal co- 
ordinate, ^ is a damping parameter, 
and K(x,x) is the ADS. 

In the initial outer-space application, 
the ADS was represented by the follow- 
ing nonlinear function: 

K(x,x) = a + bA(x,x)+ cA(x,x ) , 

where a , b , and c are constant parame- 
ters to be obtained by fitting the model 
to empirical amplitude-versus-fre- 

quency data, and A(x,x) is a modal am- 
plitude. The amplitude-versus-fre- 
quency data are obtained by means of a 
moving-window estimation technique 
in which one analyzes the instanta- 
neous vibration waveform during a 
time window of about 90 percent of the 


time-average vibration period. The am- 
plitude and frequency are taken to be 
those of a sinusoid that makes the least- 
squares best fit to the instantaneous am- 
plitude during the window (see figure). 
The window is then moved by about 2 
percent of the average period and an- 
other best-fit sinusoid is found. This 
process is repeated until a suitably repre- 
sentative sample of the vibration wave- 
form has been acquired. 

The modal amplitude is given by 


A(x,x) 


x 2 + 


K(x,x ) 


where K(x,x) is any reasonable approx- 



A Nonlinear Decaying Waveform is approximated 
with a best-fit sinusoid during a moving window. 
The resulting sinusoidal amplitude and frequency 
data are collected from fits for the entire sequence 
of window positions and used to characterize the 
frequency versus amplitude of the nonlinear wave- 
form. These frequency versus amplitude data are 
then fit to an amplitude-dependent stiffness (ADS) 
representation. 


imation of K(x,x) . One can refine the 
approximation iteratively, starting from 
K(x,x) = a, then using the resulting 
value of A(x,x) in computing a value of 
K(x,x) by use of the above equation for 
K(x,x) . 

The second model, denoted the mo- 
ment-expansion (ME) model, is embod- 
ied in the equation 

x + M(x,x) = 0, 

where the function M(x,x) is a moment 
expansion that captures damping and 
stiffness effects. The moment expansion 
is given by 

3 3 

m(x,x) = yy yy 

j = 0 i=0 

where both i and j range from 0 to 3, 
except that there is no (if) = (0,0) 
term. In the original outer-space appli- 
cation, the parameters are obtained 
from (1) modal position and velocity 
estimates obtained from Kalman-filter 
states and (2) derived accelerations. 

In a test relevant to the original 
outer-space application, the ADS and 
ME models were compared with each 
other, with a linear model, and with a 
prior nonlinear model known as the 
Duffing model. The ADS model was 
found to yield the least error. 

This work was done by Paul Brugarolas, 
David Bayard, John Spanos, and William 
Breckenridge of Caltech for NASA’s Jet 
Propulsion Laboratory. For further informa- 
tion, contact iaoffice@jpl.nasa.gov. NPO- 
41360 


@ Simpler Adaptive Selection of Golomb Power-of-Two Codes 

The selected code-parameter value is within 1 of the 
optimum value. 

NASA’s Jet Propulsion Laboratory, Pasadena, California 


An alternative method of adaptive se- 
lection of Golomb power-of-two (GP02) 
codes has been devised for use in effi- 
cient, lossless encoding of sequences of 
non-negative integers from discrete 
sources. The method is intended espe- 
cially for use in compression of digital 
image data. This method is somewhat 
suboptimal, but offers the advantage in 


that it involves significantly less compu- 
tation than does a prior method of adap- 
tive selection of optimum codes through 
“brute force” application of all code op- 
tions to every block of samples. 

A rather lengthy discussion of back- 
ground is necessary to give meaning to a 
brief summary of this innovation. For 
positive integer, m , the mth Golomb 


code defines a reversible, prefix-free 
mapping of non-negative integers to 
variable-length binary code words. 
Golomb codes are optimum for geomet- 
rically distributed sources (a model that 
frequently arises in image compression) : 
In the case of a geometrically distributed 
random variable, 5, the appropriately se- 
lected Golomb code minimizes the ex- 
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pected code-word length over all possi- 
ble lossless binary codes for 5. 

In a GP02 code, m = 2 k , where k is a 
non-negative integer. Such a code makes 
the coding process particularly simple: 
The code word for the integer 8 consists 
of the unary representation of L8/2*l (that 
is, LS/2^J zeros followed by a one) concate- 
nated with the k least significant bits of 
the binary representation of 8. More 
specifically, the code is called a GP02 
code of parameter k. 

The problem is to calculate or esti- 
mate the value of code parameter k that 
minimizes the expected bit rate (the av- 
erage number of encoded bits per 
source symbol) for an image or other 
source. This problem arises in Rice cod- 
ing, which is a coding method well 
known among experts in data compres- 
sion. The Rice algorithm encodes a 
block of samples by use of the best code 
option for the block from among several 
candidate codes that consist mostly of 
different GP02 codes. A fixed number 
of bits are used preceding the encoded 
block to indicate which code was se- 
lected. The Rice method does not spec- 



Upper and Lower Bounds on the optimum value 
of the code parameter make it possible to re- 
duce the number of code options that must be 
considered. 

ify how to find the best code option, and 
the most common approach is to ex- 
haustively try every code option to pick 
the best one for each block. Information 
from previously coded blocks is not uti- 
lized. This concludes the background in- 
formation. 

In the present method, unlike in the 
Rice method, one utilizes the mean 
sample value in each block. The 
method is based partly on a theoretical 


derivation of bounds on the optimum 
value of k as functions of the mean 
sample value (see figure). These 
bounds are such that no more than 
three code choices can be optimum 
for a given mean sample value. For a 
given mean value, one of the three 
candidate codes is selected in a proce- 
dure that involves only integer arith- 
metic (without divisions) and table 
look-ups. It has been shown that the 
value of k selected in this relatively sim- 
ple procedure is always within 1 of the 
optimum k value for the source, and 
that the cost added by the suboptimal- 
ity of the selection is never more than 
1/2 bit per sample and no more than 
about 13-percent inefficiency. In prac- 
tical image compression experiments, 
the cost added by the suboptimality of 
the selection is negligible. 

This work was done by Aaron Kiely of Cal- 
tech for NASA’s Jet Propulsion Laboratory. 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-41336. 
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Books & Reports 


® VCO PLL Frequency Syn- 
thesizers for Spacecraft 
Transponders 

Two documents discuss a breadboard 
version of advanced transponders that, 
when fully developed, would be installed 
on future spacecraft to fly in deep space. 
These transponders will be required to 
be capable of operation on any deep- 
space-communications uplink frequency 
channel between 7,145 and 7,235 MHz, 
and any downlink frequency channel be- 
tween 8,400 and 8,500 MHz. The docu- 
ment focuses on the design and opera- 
tion of frequency synthesizers for the 
receiver and transmitter. Heretofore, fre- 
quency synthesizers in deep-space 
transponders have been based on dielec- 
tric resonator oscillators (DROs), which 
do not have the wide tuning bandwidth 
necessary to tune over all channels in the 
uplink or downlink frequency bands. To 
satisfy the requirement for tuning band- 
width, the present frequency synthesizers 
are based on voltage-controlled-oscillator 
(VCO) phase-locked loops (PLLs) imple- 
mented by use of monolithic microwave 
integrated circuits (MMICs) imple- 
mented using inGaP heterojunction 
bipolar transistor (HBT) technology. 
MM1C VCO PLL frequency synthesizers 
similar to the present ones have been 
used in commercial and military applica- 
tions but, until now, have exhibited too 
much phase noise for use in deep-space 
transponders. The present frequency syn- 
thesizers contain advanced MMIC VCOs, 
which use HBT technology and have 
lower levels of flicker (1/f) phase noise. 
When these MMfC VCOs are used with 
high-speed MMIC frequency dividers, it 
becomes possible to obtain the required 
combination of frequency agility and low 
phase noise. 

This work was done by Scott Smith, 
Narayan Mysoor, James Lux, and Brian 
Cook of Caltech for NASA’s Jet Propulsion 
Laboratory. For further information, contact 
iaoffice@jpl. nasa.gov. 

NPO-42909 


® Wide Tuning Capability for 
Spacecraft Transponders 

A document presents additional in- 
formation on the means of implement- 
ing a capability for wide tuning of mi- 
crowave receiver and transmitter 


frequencies in the development re- 
ported in the immediately preceding 
article, “VCO PLL Frequency Synthesiz- 
ers for Spacecraft Transponders” (NPO- 
42909). The reference frequency for a 
PLL-based frequency synthesizer is de- 
rived from a numerically controlled os- 
cillator (NCO) implemented in digital 
logic, such that almost any reference 
frequency can be derived from a fixed 
crystal reference oscillator with micro- 
hertz precision. The frequency of the 
NCO is adjusted to track the received 
signal, then used to create another 
NCO frequency used to synthesize the 
transmitted signal coherent with, and at 
a specified frequency ratio to, the re- 
ceived signal. The frequencies can be 
changed, even during operation, 
through suitable digital programming. 

The NCOs and the related tracking 
loops and coherent turnaround logic 
are implemented in a field-programma- 
ble gate array (FPGA) . The interface be- 
tween the analog microwave receiver 
and transmitter circuits and the FPGA 
includes analog-to-digital and digital-to- 
analog converters, the sampling rates of 
which are chosen to minimize spurious 
signals and otherwise optimize perform- 
ance. Several mixers and filters are used 
to properly route various signals. 

This work was done by James Lux, 
Narayan Mysoor, Biren Shah, Brian Cook, 
and Scott Smith of Caltech for NASA’s Jet 
Propulsion Laboratory. For further informa- 
tion, contact iaoffice@jpl.nasa.gov. 
NPO-43000 


4$ Adaptive Deadband Syn- 
chronization for a Space- 
craft Formation 

A paper discusses general problems in 
estimation and control of the states (posi- 
tions, attitudes, and velocities) of space- 
craft flying in formation, then addresses 
the particular formation-flying-control 
problem of synchronization of dead- 
bands. The paper presents a deadband- 
synchronization algorithm for the case in 
which the spacecraft are equipped with 
pulse-width-modulated thrusters for 
maintaining their required states. The al- 
gorithm synchronizes thruster-firing 
times across all six degrees of freedom of 
all the spacecraft. The algorithm is scala- 
ble, inherently adapts to disturbances, 
and does not require knowledge of space- 


craft masses and disturbance forces. In 
this algorithm, one degree of freedom of 
one spacecraft is designated the leader, 
and all other degrees of freedom of all 
spacecraft as followers. The Cassini adap- 
tive optimum deadband drift controller is 
the subalgorithm for control in each de- 
gree of freedom, and the adaptation is 
run until each spacecraft achieves a speci- 
fied drift period. The adaptation is critical 
because a different disturbance affects 
each different degree of freedom. Then 
the leader communicates its thruster-fir- 
ing starting times to the followers. Then, 
for each follower, a deadband-synchro- 
nization subalgorithm determines the 
shift needed to synchronize its drift pe- 
riod with that of the leader. 

This work was done by Daniel Scharf Fred 
Hadaegh, and Bryan Kang of Caltech for 
NASA’s Jet Propulsion Laboratory. 

The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-43258. 


Analysis of Performance of 
Stereoscopic-Vision Soft- 
ware 

A team of JPL researchers has ana- 
lyzed stereoscopic vision software and 
produced a document describing its per- 
formance. This software is of the type 
used in maneuvering exploratory ro- 
botic vehicles on Martian terrain. The 
software in question utilizes correlations 
between portions of the images 
recorded by two electronic cameras to 
compute stereoscopic disparities, which, 
in conjunction with camera models, are 
used in computing distances to terrain 
points to be included in constructing a 
three-dimensional model of the terrain. 
The analysis included effects of correla- 
tion-window size, a pyramidal image 
down-sampling scheme, vertical mis- 
alignment, focus, maximum disparity, 
stereo baseline, and range ripples. Con- 
tributions of sub-pixel interpolation, ver- 
tical misalignment, and foreshortening 
to stereo correlation error were exam- 
ined theoretically and experimentally. It 
was found that camera-calibration inac- 
curacy contributes to both down-range 
and cross-range error but stereo correla- 
tion error affects only the down-range 
error. Experimental data for quantifying 
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the stereo disparity error were obtained 
by use of reflective metrological targets 
taped to corners of bricks placed at 
known positions relative to the cameras. 
For the particular l,024-by-768-pixel 
cameras of the system analyzed, the stan- 
dard deviation of the down-range dispar- 
ity error was found to be 0.32 pixel. 

This work was done by Won Kim , Adnan 
Ansar, Robert Steele, and Robert Steinke of 
Caltech for NASA’s Jet Propulsion Labora- 
tory. For further information, contact iaof- 
fice@jpl. nasa.gov 
NPO-42487 


Q Estimating the Inertia Ma- 
trix of a Spacecraft 

A paper presents a method of utilizing 
some flight data, aboard a spacecraft 
that includes reaction wheels for atti- 
tude control, to estimate the inertia ma- 
trix of the spacecraft. The required data 
are digitized samples of (1) the space- 
craft attitude in an inertial reference 
frame as measured, for example, by use 
of a star tracker and (2) speeds of rota- 
tion of the reaction wheels, the mo- 
ments of inertia of which are deemed to 
be known. 

Starting from the classical equations 
for conservation of angular momentum 
of a rigid body, the inertia-matrix-estima- 
tion problem is formulated as a con- 
strained least-squares minimization 
problem with explicit bounds on the in- 
ertia matrix incorporated as linear ma- 
trix inequalities. The explicit bounds re- 
flect physical bounds on the inertia 
matrix and reduce the volume of data 
that must be processed to obtain a solu- 
tion. The resulting minimization prob- 
lem is a semidefinite optimization prob- 
lem that can be solved efficiently, with 
guaranteed convergence to the global 
optimum, by use of readily available al- 
gorithms. In a test case involving a 
model attitude platform rotating on an 
air bearing, it is shown that, relative to a 
prior method, the present method pro- 
duces better estimates from few data. 

This work was done by Behget A$ikme§e, 
Jason Keim, and Joel Shields of Caltech for 
NASA’s Jet Propulsion Laboratory. 


The software used in this innovation is 
available for commercial licensing. Please 
contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-43631. 


© Spatial Coverage Planning 
for Exploration Robots 

A report discusses an algorithm for 
an onboard planning and execution 
technology to support the exploration 
and characterization of geological fea- 
tures by autonomous rovers. A rover 
that is capable of deciding which ob- 
servations are more important relieves 
the engineering team from much of 
the burden of attempting to make ac- 
curate predictions of what the avail- 
able rover resources will be in the fu- 
ture. Instead, the science and 
engineering teams can uplink a set of 
observation requests that may poten- 
tially oversubscribe resources and let 
the rover use observation priorities 
and its current assessment of available 
resources to make decisions about 
which observations to perform and 
when to perform them. 

The algorithm gives the rover the abil- 
ity to model spatial coverage quality based 
on data from different scientific instru- 
ments, to assess the impact of terrain on 
coverage quality, to incorporate user-de- 
fined priorities among subregions of the 
terrain to be covered, and to update cov- 
erage quality rankings of observations 
when terrain knowledge changes. When 
the rover is exploring large geographical 
features such as craters, channels, or 
boundaries between two different re- 
gions, an important factor in assessing the 
quality of a mission plan is how the set of 
chosen observations spatially cover the 
area of interest. The algorithm allows the 
rover to evaluate which observation to 
perform and to what extent the candidate 
observation will increase the spatial cover- 
age of the plan. 

This work was done by Daniel Gaines, 
Tara Estlin, and Caroline Chouinard of Cal- 
tech for NASA’s Jet Propulsion Laboratory. 

The software used in this innovation is 
available for commercial licensing. Please 


contact Karina Edmonds of the California 
Institute of Technology at (626) 395-2322. 
Refer to NPO-44282. 


$ Increasing the Life of a 
Xenon-Ion Spacecraft 
Thruster 

A short document summarizes the re- 
design of a xenon-ion spacecraft 
thruster to increase its operational life- 
time beyond a limit heretofore imposed 
by nonuniform ion-impact erosion of an 
accelerator electrode grid. A peak in the 
ion current density on the centerline of 
the thruster causes increased erosion in 
the center of the grid. The ion-current 
density in the NSTAR thruster that was 
the subject of this investigation was char- 
acterized by peak-to-average ratio of 2:1 
and a peak-to-edge ratio of greater than 
10:1. The redesign was directed toward 
distributing the same beam current 
more evenly over the entire grid andin- 
volved several modifications of the mag- 
netic-field topography in the thruster to 
obtain more nearly uniform ionization. 
The net result of the redesign was to re- 
duce the peak ion current density by 
nearly a factor of two, thereby halving 
the peak erosion rate and doubling the 
life of the thruster. (Note: NSTAR stands 
for NASA SEP Technology Application 
Readiness; SEP stands for solar electric 
propulsion.) 

This work was done by Dan Goebel, James 
Polk, Anita Sengupta, and Richard Wirz of 
Caltech for NASA’s Jet Propulsion Labora- 
tory. Further information is contained in a 
TSP (see page l).In accordance with Public 
Law 96-51 7, the contractor has elected to re- 
tain title to this invention. Inquiries concern- 
ing rights for its commercial use should be ad- 
dressed to: 

Innovative Technology Assets Management 

JPL 

Mail Stop 202-233 
4800 Oak Grove Drive 
Pasadena, CA 91109-8099 
(818) 354-2240 
E-mail: iaoffice@jpl.nasa.gov 
Refer to NPO-43495, volume and number 
of this NASA Tech Briefs issue, and the 
page number. 
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